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ABSTRACT 

The reduction of circuit equations to normal form for numerical 
integration is considered for the general circuit where excess reactive 
elements, all types of linearly dependent sources, and nonlinear dissipa- 
tive and reactive elements are present. For the linear circuit, necessary 
and sufficient conditions for the existence of numerical solutions are 
considered and stated. For the nonlinear circuit, reduction to normal 
form is not always possible. Numerical solution is shown to be simpli- 
fied if certain a priori conditions are satisfied in formulating the original 
circuit equations. A new systematic reduction procedure is presented 
for obtaining the normal form equations. This procedure is also extended 
to a new procedure for obtaining transfer and immittance functions of the 
linear circuit from a proper tree formulation . 
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I. INTRODUCTION 



The formulation of circuit equations for their numerical solution is of 

current interest. The initial contribution was made by Bashkow, who 

12 3 

formulated the concept of the proper tree. ' ' This was followed by 
the work of Bryant, and Kuh and Rohrer who refined the techniques. The 
object of their procedure, which was based primarily upon topology, 
is to obtain the circuit equations in matrix form so that one has a first- 
order differential equation. The approach uses the state-variable tech- 
nique, which has dominated the literature. However, it should be noted 
that other formulations of circuit equations, such as the TRAC Program, 
do not use state variables. 

Once the circuit equations have been formulated by the proper-tree 
concept, they have the following form 

Ax + Bw + Cz = Dx + Ew + Fz + Gu 
where x is the vector of voltages and currents associated with inde- 
pendent energy-storing elements, w is the vector of currents and voltages 
associated with dissipative elements, z is the vector of currents and 
voltages associated with excess reactive elements, u is the vector of 
independent sources, and the capital letters are coefficient matrices. 

For numerical solution it is necessary to alter the form to the normal 
form*, x = Ax + Bu for the linear case and x = f(x,u,t) for the non- 
linear case, so that numerical integration can be applied directly. 

It has been shown that the circuit matrix equation, based on the 
formulation of the proper tree, gives the relation of the currents and 
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voltages of each circuit element to the currents and voltages of all other 
circuit elements. Of this set of circuit variables, currents and voltages 
of the circuit elements, there is a minimal set called states that are 
required to describe the circuit. If those variables which do not belong 
to the minimal set can be eliminated the resulting matrix equation will be 
of normal form, or an equation that can be readily converted to normal 
form . 

Chapter II is a study of the circuit that contains linear circuit 
elements and linearly dependent sources. The formation of the proper 
tree is discussed. A new systematic technique for matrix reduction 
that eliminates the unwanted circuit variables is developed. The circuit 
that contains no dependent sources has a matrix equation that always 
reduces to normal form and reveals that the number of states is equal 
to the number of independent energy-storing elements in the circuit. 

The introduction of dependent sources alters the coefficient matrices 
of the circuit equation. These alterations make the reduction process 
conditional and these conditions are discussed. In addition, it is found 
that the introduction of dependent sources may alter the number of states 
required to describe the circuit. Those reactive variables that are not 
destined to become states are designated as surplus. It is shown that 
the surplus reactive variables can always be eliminated from the 
equations . 

In Chapter III the reduction techniques developed in Chapter II are 
utilized in establishing a new procedure for deriving driving-point 
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impedance/admittance functions and transfer functions of the linear cir- 
cuit. This procedure enables a transfer immittance to be calculated 
directly from the proper tree formulation. Derivatives are replaced by 
their frequency-domain equivalents . 

In Chapter IV nonlinear dissipative elements are included with the 
linear circuit elements. It is shown that if normal form state equations 
are to be derived the currents and voltages of the nonlinear dissipative 
elements must be restricted in their relationship to the voltages and 
currents of other circuit elements. This relationship may be controlled 
by a priori circuit modifications that increase the number of states. 

Chapter V extends the discussion to nonlinear reactive elements. 
Two new theorems relating to the conditions for the elimination of 
dissipative variables are stated. It is also revealed that the surplus 
reactive variables of the general nonlinear circuit must be associated 
linear elements and that these currents/voltages must not be related to 
the currents/voltages of nonlinear elements if normal form equations are 
desired . 
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II. REDUCTION OF LINEAR CIRCUIT EQUATIONS 
TO NORMAL FORM 



State equations are often used in circuit analysis. Normal form 
state equations, x = f(x,u,t), may be obtained from the matrix equation 
of the circuit's currents and voltages. 

A. FORMATION OF THE CIRCUIT MATRIX EQUATION 

In order to derive circuit equations in a consistent manner a proper 
tree is used.* Any tree for the circuit may be used to obtain the circuit 
equations. However, if other than the proper tree is used, additional 
steps of transposing the variables to the proper side of the equation are 
necessary in the reduction process. The proper tree is drawn as follows. 
All voltage sources are designated as tree branches and all current 
sources are designated as links. This precludes a loop of voltage 
sources, which necessarily must be constrained so that the loop can be 
eliminated, or a cut-set of current sources, which must be constrained 
so that the cut-set can be eliminated, as illustrated in Figure 1. All 
capacitive elements, less the minimum required to prevent the formation 
of a loop, are designated as tree branches (C); those exluded become 
links (elastance, S) . All inductive elements, less the minimum required 
to prevent the formation of a cut-set, are designated as links (L); those 
excluded become tree branches (inverse inductance, T ). Resistive 
elements which do not complete loops are designated as tree branches 
(conductance, G); the excess become links (resistance, R) . Thus the 
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order of precedence for drawing tree branches is (1) voltage sources, 

(2) capacitive elements, (3) dissipative elements, and (4) inductive 
elements . 

It is desirable to give consideration to elements whose currents 
and voltages are used as controlling terms for dependent sources. When 
a source is voltage dependent, the element across which the controlling 
voltage is developed should be a tree branch if this constraint does not 
upset the precedence for forming a proper tree. When a source is current 
controlled, the element through which the controlling current flows should 
be a link. Should these arrangements not be possible, the controlling 
element’s branch current may be written in terms of the link currents of 
elements forming a cut-set with the controlling element. Similarly, a 
link element's voltage may be written in terms of the branch voltages of 
elements forming a loop with the controlling element. Should the net- 
work be too complex to readily perform these substitutions and arrange- 
ments , this special consideration should be disregarded. Most circuits 
have very few dependent sources in comparison to other circuit elements. 
This consideration causes the dependency constants to appear in the 
circuit equations with the constants dictated by the network topology. 

Kirchhoff's voltage law is written for each fundamental loop defined 
by links L, R, and S. Kirchhoff's current law is written for each funda- 
mental cut-set defined by tree branches C, G, and r. All voltage sources 
appear in loop equations and all current sources appear in cut-set 
equations . Figures 2 and 3 illustrate the conceivable loops and cut-sets . 
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A loop defined by an L may contain dependent voltage sources, indepen- 



dent voltage sources , C’s , G's and T 's . A loop defined by an R may 

contain all of the foregoing with the exception of r's. A loop defined by 

an S does not contain G's and r's. A cut-set defined by a C may contain 

L's, R's, S's, and current sources. Cut-sets defined by G's do not 

contain S's and cut-sets defined by r's do not contain S’s and R's. By 

using a proper tree other loops and cut-sets are prohibited. 

In the equations which follow, F (with appropriate subscripts)'*' 

relates link and tree branch voltages and link and tree branch voltages 

and link and tree branch currents according to the topology of the circuit. 

A prime designates the transpose of a matrix, e , e , and e represent 

L K b 

independent voltage sources, j , j , and j represent independent 

r G C 

current sources, r and p (with appropriate subscripts) are submatrices of 
the proportionality constants of linearly dependent voltage sources . h 
and g (with appropriate subscripts) are the proportionality constants of 
the linearly dependent current sources, h and M are dimensionless where- 
as r and g have dimensions of ohms and mhos respectively. In the 
matrix equation which follows, fundamental loop voltage equations and 
fundamental cut-set current equations are written alternately as indicated. 
Thus, the circuit equations become (1). 
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B. THE REDUCTION PROCESS 

The reduction process includes the elimination of the variables 
associated with the dissipative elements, the transposition of the 
variables associated with the excess reactive elements, and the elimina- 
tion of the coefficient matrix from the left side of the reduced equations . 
The order of the first two steps is unimportant. 
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Partitioning (1) as indicated and lettering 
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enables (1) to be written as (3). 
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In (3), x designates the state variables, w designates the "unwanted" 
dissipative variables, and z designates the excess reactive element 
variables. Appendix A relates the submatrices of (3) to the submatrices 
of (1). 



1 . Elimination of Dissipative Variables 

The following rearrangement of (3) shows the relation of the 
variables, w, to the other variables and independent sources. 
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To have a determinate solution for the variables, w, the rank of 
the matrix multiplying w in (4) must be equal to the number of dissipative 
elements in the circuit. This matrix has dimensions n by m . Where n is 
the total number of circuit elements and m is the total number of 
dissipative elements in the circuit. Because of the dependent sources. 
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whose proportionality constants may be any real number, there exists the 
possibility that the matrix multiplying w in (4) will have rank less than 
m. In this event the solution for the variables, w, will be indeterminate. 

Let the rank of the matrix multiplying w in (4) be m-i because of 
row dependence. Rearranging the rows of (4) to group the rows such 
that the first tier is dependent upon the second tier leads to the following. 
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M ^2 is a square matrix of rank m-i and dimensions m-i by m-i. The 



enclosed subscripts denote the dimensions of the w vector partitions. 
Since the first tier of M is a linear combination of the second tier, there 
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must exist some matrix K such that K 
Multiplying the second tier of (4a) by -K and adding the results to the 
first tier yields 
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The second tier of (4b) yields the indeterminate solution for w. 

W (m-if M 22 (N 21 X + N 22 4 - T 21*- T 22 Z ' M 21 w (i) + e 2> (4c) 

Recalling that there were n elements in the circuit, m of which 
were dissipative, it can be seen that the first tier of (4b) contains n-m+i 
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equations relating to n-m reactive variables. There is no guarantee that 



solutions obtained from a set of n-m of these equations would be con- 
sistent with solutions obtained from a different set. Assuming that the 
state equations obtained from the n-m+i equations were consistent, the 
substitution of these states and their derivatives into (4c) does not yield 
a determinate solution for the dissipative variables. Therefore, it is 
necessary that the coefficient matrix of the w variables in (4) have rank 
equal to the number of dissipative elements in the circuit if the w 
variables are to be eliminated. 

The submatrix 

then the solution for w is found in the second tier of (4) 



P -Q has dimensions m by m . If the rank is m, 
22 22 J 
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Substituting (5) into the first and third tiers of (4) yields 
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where 
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In the event that (P -Q ) is singular, the matrix inverse used 

Li Zi Ci Ci 

in (5) would not exist. However, the n rows of (4) may be rearrange to 
provide a square matrix of rank m. This matrix may then be used in lieu 
of (P 2 2 — ^ 2 2^ * n an ecIuat * on s i m il ar to (5) yielding results similar to 
(6). Failure to achieve a suitable arrangement of the rows (4) leads to 
an indeterminate solution for the variables, w. Since these variables 
are resistor currents and voltages and must remain finite, indeterminate 
solutions are not anticipated in circuit analysis problems. 

2 . Transposition of Excess Reactive Variables 

To achieve the normal form equations from (6), z and z must be 
transposed. This transposition yields 
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The variables, z, are currents and voltages of reactive elements 
that were deemed excess in the formation of the proper tree. In a circuit 
free of dependent sources these variables would never become state 
variables . Because of the dependent sources (8) shows that this 
distinction can no longer be made. Thus, if 
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(8) becomes 

Ay - By + U (10) 

Further reduction depends upon the rank of the matrix A. If A 
is nonsingular, then (10) yields the normal form 



y = A 1 By + A _1 U (11) 

If the matrix A is singular, (11) is not valid and additional steps in the 
reduction process are required. 

3 . Elimination of Surplus Reactive Variables 

The matrix A in (10) can be singular because of linear dependence 
between the rows or between the columns or both. If the rows of (10) are 
regrouped forming one tier which consists of the independent rows of the 
matrix A and a second tier with the remaining rows and these tiers are 
partitioned to form square submatrices along the principal diagonals of 
A and B, the following results. 
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In (12) the subscripts i and d denote independence and depen- 
dence respectively. 

Since the second tier of the matrix A in (12) is a linear 
combination of the first tier, there exists some matrix K such that 
jj^di '^dd" 1= ^^ii ^id ’ Multiplying the first tier by -K and adding the 
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results to the second tier yields 
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Assuming that (B^-KB^) is nonsingular, the second tier of (13) yields 
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(14) 



Substituting (14) and the assoicated derivative equations into the first 
tier of (13) yields 



A*y. = B*y. + C*U. + D*U, + E*U. + F*U 



(15) 
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If A* is singular because of row dependence the process used to 
derive (15) from (10) is repeated. Each repetition of this process may 
yield higher-order derivatives of the source terms. If A* is singular 
because of column dependence a process similar to that used to derive 
(15) is employed. 

Assume that the matrix A in (10) is singular because of depen- 
dence between the columns. Regrouping and partitioning (10) 
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accordingly yield an equation identical to (12) where 
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Multiplying the first tier of (12), partitioned according to column 
dependence, by A yields 
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Multiplying "the first tier of (18) by -A^ and adding the results to the 



second tier yields (19). 
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Substituting the solution for y from the second tier of (19) 



and the associated derivative equations into the first tier yields 
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If the matrix *A is singular and has row dependence, the 
process beginning with (12) is repeated. If there is only column depen- 
dence, the latter part of the process is repeated. 

In the foregoing it was assumed that the matrices (B^-KB ^), 

used in (13) through (16), and ^d^di^i/^id^ ' used in through 

(21), were nonsingular. If these matrices are singular, there will exist 
some rearrangement of the columns of (13) and (19) that will provide 
nonsingular matrices that can be used in place of (B^-KB ^) and 

(B , -A A./b. ,). This statement is proved by contradictions that show 
dd dx 11 id 

that the rank of the second-tier submatrix of y coefficients in (13) or 
(19) is equal to the number of variables to be eliminated. First, it is 
shown that there are no dependent rows; secondly, it is shown that the 
number of independent columns is equal to or greater than the number of 
variables to be eliminated. 

Let the following represent the arrangement of (13) or (19) in 
which the independent rows of the second tier of the y coefficient 
matrix have been separated from the dependent rows, the dependent rows 
being placed in a third tier. 



A.. 

n 


'A .71 

id 


!y.l 




B.. 

li 


B id" 


~y. 




u. 

1 


0 


0 


1 

i 

i • 


= 


B* 

di 


B d*d 


! i 

( 


+ 


u d* 


0 


0 


1t<l 




B" 
_ di 




L y d_ 




_v 



Since [b" , B" I = K j~B* B* . ~j , K being the matrix that shows 
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the linear combination of the second tier of (22) that results in the third 



tier, multiplying the second tier by -K and adding the results to the third 
tier yields 
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It can be seen that the third tier of (22a) contains only equations 
of constraint on independent sources. Since these sources cannot be 
constrained and remain independent, it must be assumed that the second 
tiers of the y coefficient matrices in (13) and (19) contain no dependent 



rows . 



Assume that the second tiers of (13) and (19) do not contain a 
sufficient number of independent columns in the y coefficient matrices 
to perform the required substitutions. Let the following represent the 
rearrangement of the second tiers of (13) or (19) which provides one 
square, nonsingular submatrix, the first tiers not being shown. 
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Let Cji be the nonsingular, square submatrix. Since the second set of 
columns of C is a linear combination of the first set, there must exist 
some matrix K such that 

Multiplying the first tier of (22b) by j yields 
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Multiplying the first tier of (22d) by and adding the results to the 



second tier yields 
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Again constraint equations on the independent sources have appeared. 
Thus, it must be assumed that there are a sufficient number of inde- 
pendent columns in the second tiers of the y coefficient matrices in (13) 
and (19) to perform the operations required to eliminate the surplus 
reactive variables. Hence, the rank of the second -tier submatrix of y 
coefficients in (13) and (19) is equal to the number of variables to be 
eliminated . 



C. CONCLUSIONS 

In a circuit without dependent sources the number of states required 
to describe the circuit is equal to the number of capacitive branches plus 
the number of inductive links in the proper tree . When the circuit 
contains dependent sources, the proper tree does not necessarily 
determine the number of states. It can be seen by comparing (7) to (1) 
that additional states maybe required if the dependent sources are 
controlled by currents or voltages of reactive elements made excess by 
the formation of the proper tree. 
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Since the dependent source proportionality constants may take on 



any value, there exists the possibility that fewer states than predicted 
by the proper tree may be used to describe the circuit. Because of 
matrix singularities some of the reactive variables become surplus. The 
removal of these surplus variables introduces derivatives of the inde- 
pendent source terms and reduces the number of states. 



D. EXAMPLES 

1 . Increase in States Due to Dependent Sources 

Consider the circuit of Figure 4 which illustrates the effect of 
dependent sources that requires additional states to describe the circuit. 
The circuit equations for Figure 4(a) are the following. 
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With the absence of dissipative elements (22) may be expressed in the 
form of (6) in which 
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; C = V = 6 
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Transposing the excess reactive variables, the equations 



become 
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(24) 



(24) is in the form of (13) and the variables V and I may 

b 1 X * 

be eliminated. With the dependent sources present as shown in Figure 



4(b). 
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and the equations in the form of (24) become 
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( 26 ) 



jj) 0 0 

If 6 is nonsingular, none of the variables can be eliminated. If 
6 is singular- with rank of one, either V or I may be eliminated. 

b 1 1 \J 

Hence it can be seen that the presence of the dependent sources of 
Figure 4(b) caused an increase in the number of states required to 
describe the circuit. 

2 . Decrease in States Due to Dependent Sources 

Consider the circuit of Figure 5 which illustrates the effect of 
dependent sources that reduces the number of states by causing matrix 
singularities . The circuit equations for Figure 5 are the following . 
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The equations of (27) are converted to the form of (3). Since 
there are no excess reactive elements P 31 = Q 31 = ^33 = Q 32 = 



P 33 Q 33 



P 13 Q 13 P 23 Q 23 °* 
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From (28) it can be seen that ( P 22 _< ^22^ iS singu ^ ar anc * can 
not be used in the solution for the dissipative variables. Noting that 



P^ 2 - Q^ 2 ^ ■‘• s nonsingular, the first two rows of (28) are used to solve for 



these variables . 
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Substituting (29) into the second tier of (28) yields 
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(30) 



(30) is in the form of (13). From the first row it can it can be 

• • 

seen that V = V. Thus substituting V for V in the second row of 
(30) yields the final result. 



G 2 LI L ’ -*1 ~ G 2 CG i v 



(31) 



£rom (31) it can be seen that V is not needed to describe the 
* o 

circuit of Figure 5. Had the proportionality constant r been any other 



value, V would have been a state. Thus it can be seen that dependent 
sources may reduce the number of states required to describe the 
circuit. 
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III. LINEAR CIRCUIT EQUATIONS 
REDUCED TO IMMITTANCE FORMS 



Linear circuit equations derived by the proper-tree concept can be 
reduced to impedance, admittance, and transfer functions if the circuit 
equations are augmented by output equations . When dealing with 
transfer/immittance functions the Laplace operator with zero initial con- 
ditions must be used in lieu of derivative notation. If the circuit equa- 
tions are x = Ax + Bu and the output equations are y = Cx + Du, then x 
may be found from the circuit equations and substituted into the output 
equations to yield Y(s) = [C(sl-A) + D] U(s). This is an extended 
determination of a transfer immittance which reduces to the usual one 
when there is a single input and a single output. This procedure requires 
one to obtain the state equations and output equations and from them 
derive equations that yield the desired functions . The procedure develop- 
ed here enables one to derive the desired functions directly from (1) by 
reduction techniques. 

A. DERIVATION OF DRIVING POINT ADMITTANCE/IMPEDANCE 

If a network is driven at one set of terminals by a voltage source 
and a proper tree is drawn of the network, only one fundamental cut-set 
will be described by the source . The sum of the currents of the link 
elements forming the cut-set with the source is the current passing 
through the source. Thus, the circuit equations may be augmented with 
an output equation which expresses the source current in terms of link 
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currents . Solving for the source current in terms of the source voltage 
results in an equation of the following form . 

I (s ) = Y(s)V (s) (32) 

I(s) is the source current, Y(s) is by definition the input admittance 
function, and V(s) is the source voltage function. If V(s) is unity (32) 
becomes 

I(s) = Y(s) (33) 

If the network is driven by a unit current source, only one funda- 
mental loop will be defined by the source when the proper tree is drawn. 
The voltage developed by the source is the sum of the voltages of the 
branch elements forming the loop with the source. Thus, the circuit 
equations may be augmented with an output equation which expresses the 
source voltage in terms of the branch voltages. Solving for the source 
voltage in terms of the unit current source results in an equation of the 
following form . 

V(s) = Z(s) (34) 

V(s) is the source voltage function, and Z(s) is the input impedance 
function . 

Writing the circuit equations in the form of (1), with LI replaced 
by sLI^, etc., and performing the substitutions given by (2) yields an 
equation similar to (3), as follows. 
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Augmenting (35) with the output equation yields 
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( 36 ) 



! sP * F * F* , I 

!_ OX OW OZ J L. J 



where 0 is the output current/voltage to be determined and F , F , 

ox ow 

and F are terms which relate the output to the circuit voltages and 
oz 

currents. Th£ coupling terms, F* and F**, may include dependent 
source constants as well as topological factors. 

In the reduction process of Chapter II, to obtain the normal form 
equations, all dissipative currents and voltages were eliminated. The 
elimination process caused the dissipative parameters to appear in the 
coefficient matrices. The output equation identifies the currents and 
voltages of interest and to obtain an immittance the proper-tree equations 
should be reduced by eliminating those currents and voltages that do not 
appear in the output equation. Thus, after writing the equations in the 
form of (36), the unwanted variables are eliminated. The reduction 
process for these variables is identical to the process used to eliminate 
the dissipative variables in the derivation of normal form equations. 

(see Appendix C) 

The elimination of the unwanted variables, those for which the total 
output coupling (F** -F*) is zero, reduces the dimensions of the matrix 
equation. In addition, the reduction may cause the coefficients of the 
reduced equation to become polynomials or ratios of polynomials of the 
variables. The reduced equation is of the following form. 
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(36a) 
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where the subscript r denotes the alteration of the coefficient matrices 

and y denotes the circuit variables that were not eliminated. Any of the 

elements of P , Q , F*, and F** may now be functions of the variable s. 
r r r r 

However, this notation has been omitted for simplicity in writing. 
Solving .for y in the first tier of (36a) yields 

M = [ V ."^]' 1 [ e y j < 36b > 

Substituting this solution for y into the second tier of (36a) and solving 
for 0 yields 



0 = [>** -F*] [p -Q I 1 e (37) 

Lr rJLr rj y 

If 0 is the current of a unit voltage source, then the operations 
indicated by (37) yield an admittance function. If 0 is the voltage 
developed by a unit current source, (37) yields an impedance function. 



B. DERIVATION OF TRANSFER FUNCTIONS 

Network transfer functions may also be derived from the foregoing 
equations by letting 0 become a voltage or current at terminals not 
directly connected to the source. Care must be exercized when the out- 
put is to be an open-circuit voltage as in the case of a transfer impedance 
function. There may be elements of the circuit that are in series with the 
open-circuit terminals that do not affect the open-circuit voltage. These 
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elements should be excluded from the proper tree and the resulting 
circuit equations. 

If the network is excited by a unit current source and 0 is the open- 
circuit voltage at another set of terminals, then (37) yields a transfer 
impedance. If 0 is a current, the right side of (37) becomes a current 
gain function. If the network is excited by a unit voltage source and 0 
is a current through another set of terminals, the right side of (37) 
becomes a transfer admittance function. If 0 is a voltage, a voltage 
gain function is derived. 

By letting 0 become an output vector, both input and transfer func- 
tions may be derived simultaneously. Thus (37) may be extended to 
multiple-port circuits with multiple excitations by identifying the sources. 
Since the output equations determine the unwanted variables , it may be 
necessary to perform the elimination process in stages to obtain various 
input and transfer functions. 

C. SUMMARY 

The steps to the derivation of driving-point and transfer functions 
are summarized as follows: 

1 . determine the sources needed to derive the desired functions 

2. eliminate the variables that appear in none of the output 
equations 

3. perform the operations indicated by (37) on those output equa- 
tions where the output is coupled to all the remaining variables 
including those variables that are not coupled but not effective. 
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D. EXAMPLES 



Consider the network shown in Figure 6a. Assume that the input 



admittance, Y^, for port 1-1' is desired. Excite the circuit with a unit 
voltage source at port 1-1'. The proper tree is shown in Figure 6b. The 
augmented circuit equations are 
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(38) 



The last equation is the output equation which shows the funda- 
mental cut-set relation of the source current and I . The remaining 
circuit variables are to be eliminated. Elimination of the dissipative 
variables from (38) yields 
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Elimination of the capacitive variables from (39) yields 
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Eliminating I and performing the operations indicated by (37) 



yields 
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where 



a = 


L 2 C 1 C 2 


(41a) 
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G 1 L 2 C 2 + G 3 L 2 C 1 


(41b) 
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(4 If) 



g G l G 3 L l L 2 + G 1 R 2 L 2 C 2 + r 2 G 3 L 2 C i + L 1 C 1 + L 2 C 1 

+ L 2 C 2 (4 lg) 

h = G 1 R 2 G 3 L 2 +G 1 L 1 +G 1 L 2 +R 2 C 1 +G 3 L 2 (41h) 

i - G^ + 1 (4 li) 
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Assume that the two-port J) parameter matrix for the network of 
Figure 6a is desired. The proper tree with both ports being excited by 
current sources is shown in Figure 6c. The augmented circuit equations 



are 
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From the output equations it can be seen that I I . , and I 

LZ R1 



R3 



should be eliminated. The result of their elimination and the relocation 



of sL^ is 
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Performing the steps necessary to obtain the form of (37) yields 
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Completion of the multiplication indicated by (44) reveals the desired 
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parameters. .The results are 
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+ L^s 



(46) 



The solutions have been left in a form that allows one to check their 
accuracy by inspection. 

The hybrid parameters could be derived by exciting the circuit at 
one port with a voltage source and a current source at the other port. 
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IV. REDUCTION TO NORMAL FORM 
WITH NONLINEAR DISSIPATORS 



Because of the nonlinear current-to-voltage relations of semi- 
conductors, the numerical solutions for circuit containing these devices 
may become complex if the state variable approach is used. The circuit 
can be represented by a set of normal-form equations that are augmented 
by equations which also must be solved to obtain the nonlinear currents 
and voltages that appear in the normal-form equations. The general 
form of these equations are 



V 
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V 
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Cfj(v) + Df 2 (i) 




A 2 * 


+ 


B 2 u 


H- 


£ 3 (i,v) 


(47) 


V 


+ 


B 3 u 


+ 


f 4 (i,v) 





where the f'x are nonlinear functions and v and i are the voltages and 
currents of the semiconductors. If f. and f are zero, then (47) reduces 

O “I 

to one nonlinear differential matrix equation which may be integrated 
directly. 

4 

Diodes either by themselves, or as part of the Ebers-Moll 
equivalent for the transistor may be considered in terms of their current/ 
voltage characteristics. For the solid-state diode it is possible to 
write 

i = I (e kV d -1) 
d o 



and 



1 , /d + ^o, 
v a = k ln( — — > 



- 3-7 



The exact nonlinear relation is unimportant for the present consideration 



so that the form 

i = f(v ) (48a) 

d d 

and the inverse form 

V d = ( 48b ) 

are used here. 

5 6 

Werther and Parker and Calahan consider the diode to be a 
current-dependent voltage source or a voltage-dependent current source 
where the controlling voltage or current is written in terms of the voltages 
and currents of other circuit elements as in (47). This treatment requires 
the solution of network equations for the determination of the controlling 
voltages and currents. If the both the current function, (48a), and the 
voltage function, (48b), are single valued (the tunnel diode is exluded) , 
then the nonlinear dissipative element may be treated as the linear 
dissipative elements are treated when writing circuit equations. That 
is, the device may be used as either a branch element or a link element 
when forming a proper tree. 

A. TOPOLOGICAL CONSIDERATIONS 

The nonlinear elements are given precedence between capacitive 
elements and linear dissipative elements when establishing the branches 
of the proper tree. When writing a voltage equation for a loop defined by 
a nonlinear, dissipative link element, a term on the left side of the equa- 
tion such as RI is replaced by the voltage v . The term corresponding to 
R d 

- 3 £ 



I of the right side of the equation is f(v ,). Similarly, when writing a 
R d 

current equation of a cut-set defined by a nonlinear, dissipative element, 

the term on the left side of the equation is i , and the term on the right side 

d 

is g(i .). The matrix equation for the circuit will be similar to that of (1) . 
d 

B. ELIMINATION OF DISSIPATIVE VARIABLES 

If state equations are to be obtained the dissipative variables must 
be eliminated. While the linear dissipative variables may readily be 
eliminated, the elimination of the nonlinear variables depends upon the 
relationship of the currents/voltages of the nonlinear elements to the 
voltages and currents of other circuit elements. This may be seen by 
the following discussion. 

1 . Nonlinear Variables As General Functions of Other Variables 



Assume that the nonlinear voltages/currents are functions of all 
other voltages and currents. After writing the circuit equations and 
completing the substitutions given by (2) the equations become 
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(49) 

where x, z , and w are given by (2) and P and Q. ^ (i , j , = 1,2,3) are 
given by Appendix A. P, ^ and Q. . (i = 4,5; j = 1,2, 3, 4, 5 and J = 4,5,* 
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i = 1,2 ,3 ,4,5) are matrices which relate the voltages of nonlinear 
dissipative links to other voltages in the fundamental loops which are 
defined by nonlinear elements, and the currents of nonlinear branches to 
the currents of other elements which are in the cut-sets that are defined 
by nonlinear, dissipative branch elements. 

The first two tiers and corresponding columns of (49) may be 
combined by transposing z and z . Thus, by forming y from x and z as 
was done in (9), (49) becomes the following. 
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(50) 



As was discussed in Chapter II, the w variables may be eliminat- 
ed if the rank of the combined coefficient matrices for these variables is 
equal to the number of these variables. Their elimination reduces (50) 
to the following . 
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This 



The multiplication of (51) by [K] ^ yields (47). Further 
reduction depends upon the elements of M ^ ^ 32 ' anc * ^ 
may be seen by writing the second two tiers of (51) as follows. 



V d K 22 K 23 



*d K 32 K 33 



-1 



91 *'*o, V o M M f(v ) 

21 [y] + 21 [y] + 2 + 22 23 d (53) 



31 
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31 



V 3 M 32 M 33 g(i d } 



Equation (5 3) suggest several undesirable solutions for the 
nonlinear currents and voltages . Transcendental solutions may be 
required. If .these are not required, solutions which reveal that a non- 
linear current/voltage is a nonlinear function of other currents/voltages 
may result. The substitution of the latter solutions into the first tier 
of (51) yields functions of functions. If M22 = ^23 = ^32 = ^33 = 13 ' 



these undesirable solutions cannot exist. By comparing (52) with (50), 



it can be seen that the desired conditions prevail if ^23 = ^24 = ^32 = ^33 = ^34 



= 0. This implies that the fundamental loops defined 
by nonlinear elements must not contain other dissipative elements, 
either linear or nonlinear. Similarly, cut-sets defined by nonlinear 
elements must not contain other dissipative elements. These restric- 
tions must also apply to dependent sources which are controlled by the 
currents/voltages of dissipative elements. Thus, if desirable equations 
are to result, the nonlinear, dissipative variables must be restricted to 
being functions of reactive variables only. 



2 . Nonlinear Variables As Functions of Reactive Variables 

If the currents and voltages of nonlinear elements are not allowed 
to be controlling terms for dependent sources and these currents/voltages 
can only be related to capacitor voltages or inductor currents , then (50) 
becomes the following. 
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The eliminations of the w variables yields the following. 

xx x, 

0 
0 



K n 0 0 I 



I 0 I 

0 I 



v di = 



J *d_ 



M , 


M 


11 






0 


21 




M„ , 


0 


31 




. and 


i , 


d 


d 



y 




~ v r 


f(v d ) 


+ 

1 


V 2 


5 (1 d! 


[ 

1 


: v 3. 



(55) 



K ny - M n y + W 1 + M 12 f(M 21 y + V 2 ) + M 13 g(M 31 y f Vg) 



(56) 



If in (56) is nonsingular, then (56) readily yields the normal 
form equations. If K^, is singular some of the variables in y must be 
eliminated. The elimination of the surplus reactive variables may be 
possible but very difficult. Assume that the left side of (56) is zero 
and that one of the variables y is to be written in terms of the other 
variables and sources. It can be seen that the solution of a transcen- 
dental equation may be involved. To investigate the conditions neces- 
sary for reduction when is singular; a circuit that is free of 
dependent sources is considered. 
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For a circuit that does not contain dependent sources the circuit 



equations reduce to the following. 
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and Q will be other than zero 

0 o 



whenever a fundamental cut-set is formed linear dissipative elements 

and nonlinear dissipative elements. Likewise, Q and Q ^ will be 

other than zero whenever a fundamental cut-set is formed with nonlinear 

dissipative elements. The elimination of the linear dissipative variables 

gives rise to coefficients for f(v ,) and g(i .) in the reduced equations 

d d 

causing the solutions for v . and i , to become transcendental. If i . is 

d d d 

eliminated from (57), then Q and Q cause the solution for v to 

T 0 D T Cl 

become transcendental. Thus, normal-form state equations will not 
result from the reduction of circuit equations when fundamental cut-sets 
are formed with nonlinear dissipative elements and other dissipative 
elements . 

If the matrices mentioned in the preceeding paragraph were zero, 
the nonlinear voltages/currents would be functions of the state variables 
and sources only and all dissipative variables could be eliminated. The 
resulting normal form equations would be as follows. 

x = Ax + BU + CU + Df(Ex + FU) + Gg (Hx + JU) (58) 
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C. CONCLUSIONS 



The circuit equations for circuits containing nonlinear dissipative 
elements do not readily reduce to normal-form state equations . By a 
priori circuit modifications of placing an inductor in series with those 
nonlinear elements which have loop equations formed with variables that 
are not destined to become state variables and placing a capacitor in 
parallel with those nonlinear elements which have cut-set equations 
formed with variables that are not destined to become state variables, 
the nonlinear elements' voltages and currents are forced to become state 
variables . While these modifications do give rise to additional states 
they allow the circuit equations to be reduced to normal-form equations . 



D. EXAMPLE 



Consider the circuit shown in Figure 7a . The circuit equations are 
the following . 
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(59) 



This equation cannot be put in normal form because of the coupling 
between diodes. If the circuit is modified as in Figure 7b, the circuit 
equations become 
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The diodes are now uncoupled and (60) may be placed in normal 



form . 
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(61) is obtained by eliminating from (60). 
are obtained by multiplying (61) by 

~C 0 61 -1 

0 C 0 

0 0 L 



Normal-form equations 



(62) 
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V. REDUCTION TO NORMAL FOR M 
WITH GENERAL NONLINEAR ELEMENTS 



When nonlinear inductors and capacitors are present, the circuit 
equations may be reduced to the following form. 
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(64a ,b) 



(65a, b) 



0 is the flux associated with inductors and Q is the charge 

G 

associated with capacitors . 0 and Q are functionally related to the 

J-j c 

inductor currents and capacitor voltages as given by (64b). For the 

linear case 0 T = LT and Q =CV . (65a, b) give the general relationships 

L L c c 

for excess reactive elements. 

If x and z in (49) are replaced by f (x) and f (z) respectively, the 

X z 

result is (63). Since z and f z (z) may be transposed, it is possible to 
form 
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(66a, b) 
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If f (y) is used in place of y in the equations that begin with (50), the 
arguments for reduction with nonlinear inductors and capacitors would 
be obtained. Thus, if the dissipative variables are to be eliminated, 
the conditions for the elimination of these variables from (49) must also 
apply to (63) . 

The conditions for the elimination of the dissipative variables are 
stated in the following theorems. 

Theorem 1 . If a coefficient submatrix for the linear dissipative variables 
in the equations which relate these variables to the other circuit 
variables and sources has rank equal to the number of linear dissipative 
elements, the linear dissipative variables may be eliminated and normal- 
form equations may be derived. (For proof see Chapter IIB1) 

Theorem 2 . If each of the voltages or currents of nonlinear dissipative 
elements, nonlinear dissipative variables, are linearly related to only 
independent sources and the reactive voltages and currents that are 
functions of variables that are destined to become state variables, then 
the nonlinear dissipative variables may be eliminated and normal form 
equations may be derived. (For proof see Chapter IVB2) 

Having eliminated the dissipative variables from (63), the resulting 
equations are as follows . 

Ay = Bf (y) + U + CfjtDf (y) + U) + Eg d (Hf y (y) + U) (67) 
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If A in (67) is singular there exist surplus states to be eliminated. 



Since these states do not appear in functions f and g^ which result 



from the elimination of the nonlinear dissipative variables, an equation 
of the following form must be solved to eliminate the surplus state 
functions f (y) and the surplus state derivatives. 



where the subscripts i and s denote the state functions to be retained 
and the surplus state functions respectively. While the surplus state 



sources, the solution for the surplus state derivatives involves partial 
derivatives of both the surplus functions and the functions to be retained. 



From (69) it can be seen that if the surplus states are to be elimin- 
ated the surplus state functions must be linear and must be related to 
retained functions that are also linear. Otherwise, the partial deriva- 
tives must be entered into the coefficient matrix of the retained 
derivatives and the resulting equations become nonlinear differential 
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equations that cannot be placed in normal form. 
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APPENDIX B 







(a) Voltage Source Loops 






(b) Current Source Cut— Sets 

FIG I. ELIMINATING VOLTAGE SOURCE LOOPS 
AND CURRENT SOURCE CUT-SETS 
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FIG 2. CONCEIVABLE LOOPS 



52 










> 




INDEPENDENT 

SOURCES 



Q)--gV s 

@--hS" , V s 







hi 



L 



__.0_.gr 1 I L 

___0__hl r 



FIG 3. CONCEIVABLE CUT-SETS 






(b) Circuit With Dependent 
Sources 



FI G 4. EXAMPLE I. INTRODUCTION OF STATES 
BY MUTUALLY COUPLED EXCESS 
REACTIVE ELEMENTS 
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FIG 5. EXAMPLE 2. REDUCTION OF STATES 
DUE TO DEPENDENT SOURCE 
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FIG 6. EXAMPLES 3 a 4. IMMITTANCE FORMS 
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FIG 7. EXAMPLE 5. INTRODUCTION OF ADDITIONAL 
STATE TO OBTAIN NORMAL FORM REDUCTION 
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APPENDIX C 

Elimination of Unwanted Variables 



Consider the following equation and the elimination of the variable 



w . 



A 11 A 12 A 1 3 ■ 
A 21 A 22 A 23 
A 31 A 32- A 33 



Fjx) 



w 



B 11 B 12 B 13 
B 21 B 22 B 23.: 
B 31 B 32 B 33 



X 


e 






X 


w 


+ 


e 






w 


F (z) 




e 


z 




z 



If (^22 -B 2°^ * s nonsin< 3 u l ar ' second equation may be written as 

W = ^ A 22 -B 22* ^ B 21 X + B 23 F z^ ~ A 21 F x^ _A 23 Z + °w* 

Substituting this value for w into the first and third equations and 
combining terms yields 



A 11 ( A 12 _B 12^ A 22 -B 22^ A 2l| A 13 ^ A 12 B 12^ A 22 B 22^ A 23 

A 31~^ A 32~ B 32^ A 22 -B 22^ A 2b' A 33~^ A 32” B 32^ A 22 _B 22^ A 23 

B 11 - ^ A 12 -B 12^ A 22 _B 22^ B 2l! B 13 _ ^ A 12~ B 12^ A 22 _B 22^ B 23 

B 3l”^ A 32” B 32^ A 22~ B 22^ B 2 1 1 B 33~ ^ A 32~ B 32^ A 22 _B 22^ B 23 



°x (A 12 B 12 )(A 22 _B 22 ) G w 
e z"^ A 32 _B 32^ A 22 _B 22^ G w 



F (x) 

X 



X 



Fz< 2 > 



The pattern for the multiplication and combining of terms is readily 



established . 
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